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The effect of strong quantising magnetic field on low density stellar matter is 
investigated using Thomas-Fermi-Dirac (TFD) model. The Wigner-Seitz cell struc- 
ture is assumed for the low density matter. The significant changes in the properties 
of such low density matter in presence of strong magnetic fields are discussed. It is 
seen that the decay time scale for magnetic field decreases by at least two orders of 
magnitude in such model calculation. 
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1. INTRODUCTION 



The study of low density stellar matter at the crustal region of a neutron star or in the case 
of white dwarf matter in presence of a strong magnetic field is extremely interesting both from 
the academic as well as from the Astrophysical point of view. The recent observational data of 
magnetars indicate the possibility of very strong surface magnetic field (up to 10 15 G) jl] in some 
pulsars. The observed soft gamma repeaters (SGR) discovered in BATSE j|] and KONUS (HQ 
experiments and X-ray source observed by ASCA and RXTE show strong surface magnetic field up 
to 10 15 G. These objects are called magnetars. They pose a great challenge to the existing models 
of magnetic field evolution since they require a very rapid field evolution in isolated neutron stars. 
Now the neutron star crust plays most important role in the evolution of neutron star magnetic 
field ||. The neutron star crustal matter (at relatively low density) exhibits conventional lattice 
structure, which can be approximated by a regular arrangement of Wigner-Seitz cells, with a 
positively charged nucleus at the centre surrounded by spherical distribution of electron cloud. In 
the case of white dwarfs such approximation is also valid at low density region. In the non magnetic 
case, this kind of crystalline structure have already been studied in Solid state physics / Atomic 
Physics HQ- Equation of state for low density stellar matter is also obtained using TFD model 
in the non magnetic case || . Since the shell effect of orbital electrons significantly suppresses the 
statistical effect, it is not recommended to obtain equation of state for very low density solids in 
the laboratory using TFD model (e.g. metallic iron in the laboratory). On the other hand, in the 
case of stellar matter, since the density is high enough, such an issue does not arise. It was also 
argued that application of TFD model is valid for density < 10 gm/cc. For higher densities, since 
electrons no longer remain bound within the cells, one uses Chandrasekhar's ideal gas results for 
electrons with Coulomb lattice correction 

If magnetic field of strength ~ 10 14 -10 15 G exists at the outer crust of a neutron star or inside a 
magnetised white dwarf (!) with Wigner-Seitz cell structure of the matter, the charge distribution 
of the electron cloud within such cells must be affected significantly. As a consequence, the size of 
Wigner-Seitz cells will change. This change in cell volume should affect the equation of state of 
low density stellar matter and reduces the width of outer crust of neutron stars. If such modified 
equation of state is considered for the neutron star crust, it will significantly affect the gross 
properties including the mass-radius relation of the star. To the best of our knowledge, such an 
important issue has not been discussed before. 

In this article, we assume that the matter consists of fully ionised iron nuclei and are at rest 
at the centre of Wigner-Seitz cell. Since the matter is at relatively low density, the electrons 
surrounding the positive ions are assumed to be non-relativistic. It is further assumed that in the 
microscopic scale the magnetic field is constant and is along z-axis. As we have noticed that the 
qualitative nature of the results do not change if we consider carbon or oxygen instead of iron as 
the constituent of the matter. 

The paper is organised in the following manner. In the next section we shall develop the basic 
formalism for Thomas-Fermi-Dirac model in presence of a quantising magnetic field. We have 
concluded our results and discussed the future perspective of this work at the last section. We 
have presented a brief outline of the derivation of electron-electron exchange interaction (Fock) 
term in presence of a strong magnetic field in the Appendix A. 



2. TFD MODEL IN PRESENCE OF A STRONG MAGNETIC FIELD 

In presence of a quantising magnetic field of strength B, the electron number density is given 

by §, 

n e = —Pf (1) 

where pp is the electron Fermi momentum, and e is the magnitude of electronic charge. In the 
Thomas-Fermi model for statistical treatment of atomic structure, it is assumed that within the 



2 



Wigner-Seitz cell, the electrons move in a slowly varying spherically symmetric potential V(r). 
Then the Fermi energy fj, of an electron is given by, 

M = -eV(r) + j| (2) 

where m is the electron mass. The Fermi energy /i is independent of r, otherwise electron would 
migrate to a region of smaller p. In TFD model the electron Fermi energy is given by, 

2 

/i = ecf> — u ex (pp) — constant (3) 

where u ex is the exchange part of electron-electron interaction. 
In the non magnetic case 

e 2 

U ex {PF) = —rPF (4) 

whereas in the case of a quantising magnetic field, if all the electrons are assumed to be at the 
lowest Landau level, the exchange energy is given by, (see Appendix A) 

u ex (p F ) = a(l - e-^) (5) 

where the parameters a and (5 are functions of magnetic field strength and are given in Table I. 
Rearranging eqn.(3) in the form (see also 

|£L + ae ~0VF =fl * +e ^ ( 6 ) 
2m 

where ji* = p + a is the modified form of Fermi energy of the electron. From eqn.(6) one can 
express Fermi momentum pp as a function of p* + e<f>. The numerically fitted functional form is 
given by a simple power law, 

PF = C(n* + e<f>y (7) 

where, C and 7 are constant parameters for a given magnetic field strength. In Table I, we have 
shown the variation of C and 7 with the magnetic field strength B. The potential cj> (which is the 
direct interaction term between electron-nucleus and electron-electron) is given by the Poission's 
equation (|], 

V 2 </> = 47ren e + nuclear contribution (8) 

Since the nuclear contribution is a delta function about the origin, we can omit it for r > and 
impose the boundary condition, 

lim r4>{r) = Ze (9) 

Where Z is the atomic number (= 26 for iron). The boundary condition at the cell wall of radius r s 
is that the electric field vanishes (neutral cell condition), which gives in spherical polar coordinate, 

f r \r=r a =0 (10) 

Now, using the empirically fitted form of pp given by eqn.(7), and writing the radial coordinate r 
in the scaling from r — ax, we have the Poisson's equation (from eqn. (7)) 

g=x-V (11) 
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Where 



[i* + e<j) = —u(r) (12) 
r 



and 



*2 

irn c 



(13) 



ACBe 2 ^+ 1 Zi- 1 

The boundary condition at the cell boundary (eqn.(10)) gives, 

^ = - (14) 

ax a; 

for a; = x s . Since 7 < 1 for the whole range of magnetic field strength (10 14 G < B < 10 17 G) of 
Astrophysical interest, unlike non-magnetic case, the Poisson's equation (eqn.(ll)) does not have 
singularity at the origin. Therefore, the numerical method prescribed by Feynman, Metropolis 
and Teller |§ is not necessary in the quantising magnetic field case. This qualitative change in 
the form of Poisson's equation comes from the modified form of phase space integral of electron 
number density in presence of strong magnetic fields. The non-quantising magnetic field therefore 
can not make any qualitative change in the form of differential eqn.(ll) or in other wards in the 
electron distribution within the cell. Standard fourth-order Runge-Kutta method has been used 
to obtain numerical solution of eqn.(ll). In this case the initial value for the derivative 

are chosen by shooting method to match the boundary condition at the cell surface for three 
different magnetic field strengths, B = 10 14 G, 10 15 G, and 10 17 G. The values for x s , the surface 
scaling parameter are given in table I for the above three magnetic field strengths. As we have 
noticed, the cell radius r s = ax s decreases with the increase of magnetic field strength and are 
about an order of magnitude smaller than the non magnetic value [pi . This squeezing of Wigner- 
Seitz cell in presence of strong quantising magnetic field is analogous to the well known magneto- 
striction phenomenon observed in classical magneto-statics. The variation of u(x) with x for a 
given magnetic field strength is given by the numerically fitted functional form 

u i x ) = F77 TT (16) 

l + cxp{e(a;-x )} 

where, u$, £, xq are constant parameters for a given magnetic field strength. The variation of 
these parameters with magnetic field strength are shown in Table I. In presence of strong quatising 
magnetic field, the variation of u with x is entirely different from the non magnetic case. The 
variation is more or less like the radial distribution of matter in neutron stars. Now the equation 
of state of such cold degenerate low density matter is that due to nucleons and electrons present 
in the system. The pressure contribution mainly comes from the electrons. The nuclei are at 
rest at the centre of each Wigner-Seitz cell, therefore we can ignore their contribution to kinetic 
pressure of the system. On the other hand, the energy density of the system is mainly dominated 
by rest mass of the ions. The energy contribution from the electronic sector is about five-six orders 
of magnitude less than the ionic part, therefore one can discard the energy contribution in the 
equation of state from electronic part. The mass density p is simply given by the rest mass of 
nucleons inside the cell, then we have 

3Am B , , 

P= W^f (17) 

where tub — 1.66057 x 10 _24 g, the effective nucleon mass. In Table I we have shown the variation of 
matter density with the magnetic field strength B. Since the radius of Wigner-Seitz cells decrease 
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with the increase of B, there is an increase in mass density with the increase in magnetic field 
strength. The expression for kinetic pressure from electron sector is given by 
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p% , n . / 1 \ a 

_ + aexp(/3 PF )(p F + -j-- 



(18) 



where the Fermi momentum pp has already been expressed as a function of dimcnsionlcss surface 
parameter x s (see eqn.(7)). Eqns.(17) and (18) give the equation of state P = P(p) in terms of 
the surface parameter x s . 

In fig.l we have shown the equation of state of such low density matter in presence of strong 
magnetic fields for three different cases: the upper curve is for B = 10 14 G, middle one is for 
B = 10 15 G, and the lower one is for B = 10 17 G. As we can see from the figure that the softness 
of the matter increases with the increase of magnetic field strength. Which also means that the 
matter becomes energetically more stable with respect to non-magnetic case. We have further 
noticed from the upper and the middle curves that at high density since a large number of Landau 
levels are populated for electrons, which is equivalent to the non-quantising picture of external 
magnetic field, these two curves coincide at high density. The effect of magnetic field is completely 
washed out completely at very high density. 
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FIG. 1. Equation of state for low density neutron star crustal matter in presence of strong magnetic 
fields: B = 10 14 G (upper curve), B = 10 15 G (middle curve) and B = 10 17 G (lower curve) 



3. CONCLUSION 



The outer crust of a neutron star, particularly in the case of a strong magnetic field (magnetars ?) 
plays a crucial role in the evolution of pulsar magnetic field. It is really a great challenge to explain 
field evolution in these strongly magnetised objects using existing models of field evolution. These 
objects require a very rapid field evolution. Now the TFD model for low density matter in presence 
of strong magnetic fields shows an over all contraction of the outer crust. Since the Ohmic decay of 
magnetic field in a conducting material depends on the thickness of the region, a decrease in width 
of the outer crust by an order of magnitude will cause a rapid decay of magnetic field (at least two 
orders of magnitude decrease in decay time scale). The equation of state curves (fig.(l)) indicates 
that electrons within the Wigner-Seitz cells are more strongly bound to the positively charged 
nuclei in presence of strong quantising magnetic fields than the non-magnetic (or non-quantising) 
case. Such strong binding of electrons within the cells may decrease the electrical conductivity 
of the matter. Which will further reduce the time scale for Ohmic decay of magnetic field in the 
outer crust of these strongly magnetised stellar objects. 
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APPENDIX A: 



Exchange energy is given by (with h = c = 1) 



j=i •/ \r r 



(Al) 



For the zeroth Landau level, the wave function ip( r ) f° r electron is given by 



\J LyL z 



eB 

7T 



1/4 



exp 



2 V eB' 



exp[i(p v y + p z z)] 



(A2) 



In this particular case the sum over j can be replaced by L y L z dp' y dp' z , where L y and L z are the 
linear dimensions of the box along y and z directions respectively. The volume element d 3 r' = 
dx'dy'dz' . 

Then following Lee |L2] , we have 



dy'dz' 



1 



■exp[-i(p y -p')(y - y') - i(p a - p' z )(z - z')] 



47T 



exp {—K | a; — x' \) 



(A3) 



where IT = yj (p y - p' y ) 2 + (p z - p' z ) 2 . 

Similarly d 3 r — dxdydz. The integral J dydz = L y L z . Then we have 



1 f eB 



2 
exp 



— ] Aire 2 [ dp'dp'dxdx' exp i—K I x — x' 

7T / J y z 2K 



eB 



_ £w V 
eB) 



eB 



PjlY 

eB) 



eB 



(A4) 



To evaluate the integrals over x and x' , we change the integration variables to X and Y, where 
X = x - x' and Y = (x + x')/2. 
Now 



SZ dx exp { - k 1 ^ |} exp (-T x2 ) = exp (S) erfc r h 



\V2e~B 



(A5) 



where erfc(x) is the complimentary error function. 
Then we have 



1 



2tt 



e B I — dp dp z dY\ — exp 



K 2 
2eB 



erfc 



( K 



\V2e~B 



eB / , 
2 



exp 

The Y integral is given by 

cZFexp 



2p 



2p 12 Ap y Y <ip' y Y S 



e 2 B 2 e 2 B 2 



eB 



eB 



eB 



2Y 



Py + Py 
eB 



2eB 



(A6) 



(A7) 



Then we have after changing the integration variables from p' y and p' z to P y — p y — p' y and 

Pz =Pz~ p' z 



G 



U ex = e\j d P yd P z7 ^==e^ (|L) crfc (J^j (A8) 

where the limit of P y is from — oo to oo and P z is from to 2pp for p z = pp. 
Again putting P y = P z tan#, we have 

r 2 PF r^/ 2 / I p I \ / P 2 \ 

U ex = e 2 ir dP z sec (9 d6 erfc -L^=L sec exp N-^ (A9) 
Jo Jo \v2eS / \ 2ei V 

These double integrals have been evaluated numerically as a function of Fermi momentum . The 
fitted functional form of U ex is given by 

U ex = a[l - exp(-0PF)] (A10) 
where the parameters a and /3 vary with magnetic field strength i? and are shown in Table I. 
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Table I 



B (Gauss) 


10 i4 


10 ib 


10 iV 


a (MeV) 


0.568 


1.796 


17.909 


MeV" 1 


3.412 


1.067 


0.109 


7 


0.506 


0.527 


0.658 


C 


0.973 


0.870 


0.386 


X s 


3.096 


3.170 


4.404 


r s (A) 


0.402 


0.203 


0.123 




-0.938556 


-0.937365 


-0.936123 




1.633 


1.651 


1.944 




2.097 


2.071 


1.755 




0.213 


0.204 


0.031 


p (gm/cc) 


72.79 


572.29 


962.14 
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